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Abstract 

A comprehensive study of the color singlet heavy quark states J /if), T and {l^bb 
above T c is given. We use the Field Correlator Method (FCM) for nonperturbative 
QQ potentials and the screened Coulomb potential with the T-dependent Debye 
mass computed up to two loops in the deconned phase of QCD. We calculate binding 
energies and melting temperatures of heavy mesons and baryons in the deconfined 
phase of quark-gluon plasma and the J /if) and T disintegration cross sections via 
the gluon absorption. 
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X> ! PACS numbers: 11.25 Sq, 12.38Lg, 12.38.Mh, 14.40.Pq, 25.75.Nq 

^ : 1 Introduction 

Qh- Early applications of in-medium heavy-quark potentials have employed a phenomenolog- 
ical ansatz to implement color-screening effects into the one-gluon-exchange (Coulomb) 
potential pQ . The temperature dependence of this potential is encoded in the Debye mass 
m£)(T). If Debye screening of the Coulomb potential above the temperature of decon- 
finement T c is strong enough, then J /if} production in A+A collisions will be suppressed. 
Therefore the gold-plated signature of deconfinement in the Quark-Gluon Plasma (QGP) 
was thought to be the J /if) suppression [2J. Indeed, applying the Bargmann condition [3] 
for the screened Coulomb potential at T = T c 

m 

g: v (r) = --• ^ • e -°ev>r (1) 

we obtain the simple estimate for the number of, say, cc 5-wave bound states 

oo 

n<^J \V c (r)\rdr = ^-^f— (2) 
J 3 m D (T c ) 

o 

where fi c is the constituent mass of the c-quark and for the moment we neglect the r- 
dependence of a s . Taking \i c = 1.4 GeV and a s = 0.39, we conclude that if m D {T c ) > 0.7 
GeV, there is no J /if) bound state. Parenthetically, we note that for the potential ([T]) no 
light or strange mesons (/i ~ 300 — 500 MeV) survive . But this is not the full story. 

There is a significant change of views on physical properties and underlying dynamics 
of quark-gluon plasma (QGP), produced at RHIC, see, e.g., [4J and references therein. 
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Instead of behaving like a gas of free quasiparticles - quarks and gluons, the matter created 
in RHIC interacts much more strongly than originally expected. It is more appropriate 
to describe the nonperturbative (NP) properties of the QCD phase close to T c in terms 
of the NP part of the QCD force rather than a strongly coupled Coulomb force. 

In the QCD vacuum, the NP quark-antiquark potential is V = or. At T > T c , 
a = 0, but that does not mean that the NP potential disappears. In a recent paper [5] 
we calculated binding energies for the lowest QQ and QQQ eigenstates (Q = c, b) above 
T c using the NP QQ potential derived in the Field Correlator Method (FCM) [6] and 
the screened Coulomb potential. The Debye radii were calculated for pure gluodynamics 
using Eqs. (28), (29) of Ref. [7J and the parameters given there in. In the present paper 
we refine the results of [5] and, in particular, extend our analysis to the case of the running 
a s {r) [8]. We also calculate the dissociation cross sections of J/ip and T in collisions with 
gluons. Here as in Ref. [5] the Debye mass is evaluated in quenched QCD (the number of 
light flavors rif = 0, T c = 275 MeV), the similar results using the Debye mass evaluated 
for two light flavors (T c =165 MeV) are presented in Ref. [9]. 



2 Field correlator method as applied to finite tem- 
peratures 

The NP quark-antiquark potential can be studied through the modification of the correla- 
tor functions, which define the quadratic field correlators of the nonperturbative vaccuum 
fields: 

< trF^(x)$(x,0)F A(T (0) >= A^-xa D(x) + B^ u; xa D 1 (x), 

where A^ u -xa and B^-xa are the two covariant tensors constructed from g^ v and x^x v 
[6], $(x,0) is the Schwinger parallel transporter, x is Euclidian. At T > T c , one should 
distinguish the color electric correlators D E (x), Df(x) and color magnetic correlators 
D H (x), D^(x). Above T c , the color electric correlator D E (x) that defines the string 
tension at T = becomes zero [TU] and, correspondingly, a E = 0. The color magnetic 
correlators D H (x) and D^(x) do not produce static quark-antiquark potentials, they only 
define the spatial string tension a s = a H and the Debye mass rrid oc ^Jd~ s that grows with 
T. 

The main source of the NP static QQ potential at T > T c originates from the color- 
electric correlator function D E (x): 

l/T r 

K P (r,T) = J dv{l-vT) J \d\Df(x). (3) 



In the confinement region the function Df (x) was calculated in [TTj : 

Df(x) = B eXp( " MoX) , (4) 

x 

where B = Qa^ajMo, a{ being the freezing value of the strong coupling constant to be 
specified later, <jy is the string tension at T = 0, and the parameter M has the meaning 
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of the gluelump mass. In what follows we take <jf = 0.18 GeV 2 and M = 1 GeV. Above 
T c the analytical form of should stay unchanged at least up to T ~ 2T C . The only 
change is B — > B(T) = £(T)B, where the T-dependent factor 



M T-T c 
1-0.36-f - 
B T r 



(5) 



is determined by lattice data [12]. Integrating Eq. (j^J) over A, one obtains 



V np (r,T) 
where 

and 



B(T) 
M 



l/T 



V(oo,T) 



vT) (e 
B(T) 



-uM 



Ml 



T 

Mo 



exp 



Mo 
T 



(7) 



l/T 



V(r, T) = ^21 [ (1 - vT) exp(-v / z/ 2 + r 2 M )du. 
M J 
o 

The approximate expression 



V(r,T) 



B{T) 
Ml 



T 

Ki{x)x - — exp(-x)(l + x] 

Mq 



(9) 



where x = M$r and K\(x) is the McDonnald function, has been used in [5]. For T = T c 
expressions ©and ([9]) are almost indistinguishable. At T > T c the exact potential is 
slightly more attractive, the difference between OH]) and OH]) increases with T. Even more 
drastic approximation V(r, T) oc xK\(x) was proposed [TT]. This approximation was used 
in [5] to calculate the QQQ states. 



3 Coulomb potential 

We use the perturbative screened Coulomb potential (CQ) with the r-dependent QCD 
coupling constant a s (r,T). Note that in the entire regime of distances, for which at 
T = the heavy quark potential can be described well by QCD perturbation theory, 
a s (r,T) remains unaffected by temperature effects at least up to T < 3T C and agrees 
with the zero temperature running coupling a s (r, 0) = a s (r). For our purposes, we find 
it convenient to define the r-dependent coupling constant in terms of the q 2 - dependent 
constant a#(q 2 ) calculated in the background perturbation theory (BPTh) [13J: 



oc 



2 f l sin qr 



a s (r) = - dq a B (q )• (10) 

Tr J q 



o 



The formula for a#(q 2 ) is obtained by solving the two-loop renormalization group equation 
for the running coupling constant in QCD 

/ <a 47t / Bi lnt\ . q 2 + m 2 R 

« S (q 2 ) = -r- 1 - -^7— , t = In 4 A2 B , (11) 



0ot V A) 2 t J ' ~ A 2 , 
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where are the coefficients of the QCD /3- function. The parameter ~ 1 GeV has the 
meaning of the mass of the lowest hybrid excitation. The result can be viewed as arising 
from the interaction of a gluon with background vacuum fields. Note that as(r) increases 
with Ay and, for fixed Ay, decreases with tub- 

We employ the values Ay = 0.36 GeV, mg = 0.95 GeV, which lie within the range 
determined in [T4]. The result is consistent with the freezing of a_e( r ) with a magnitude 
0.563 (see Table 1 of [15] ). The zero temperature potential with the above choice of the 
parameters gives a fairly good description of the quarkonium spectrum. 

The Debye mass mo(T) in Eq. ([I]) is expressed in terms of the spatial string tension 
<J S (T) due to chromomagnetic confinement: m£,{T) = 2.06a/<t s (T). The latter has been 
computed nonperturbatively up to two loops in the deconfined phase of QCD [7]. As was 
stated above, in this paper, we consider the pure-gauge SU(3) theory (T c = 275 MeV), 
for which mo varies between 0.8 GeV and 1.4 GeV, when T varies between T c and 2T C . 



4 Results 

4.1 Quark-antiquark states 

In the framework of the FCM, the masses of heavy quarkonia are defined as 

m n 

Mqq = + [iq + E (m Q ,tJ,Q), (12) 

Eo(rriQ,fiQ) is an eigenvalue of the Hamiltonian 

H = H + K P + V c , (13) 

rriQ are the bare quark masses, and einbeins fii are treated as c-number variational pa- 
rameters. The eigenvalues E (rrii, m) of the Hamiltonian f[T3"j) are found as functions of 
the bare quark masses and einbeins /ij, and are finally minimized with respect to the 
//j. With such simplifying assumptions the spinless Hamiltonian Hq takes an apparently 
nonrelativistic form, with einbein fields playing the role of the constituent masses of the 
quarks. Once ttiq is fixed, the quarkonia spectrum is described. The dissociation points 
are defined as those temperature values for which the energy gap between V(oo,T) and 
E (T) disappears. 

In our calculations, we use the quark-antiquark potentials whose parameters are listed 
in Table [TJ The potential I was employed in [5] . This potential uses the approximation 
© for V(r) and the constant value a s = 0.35 for the Coulomb potential. The potential 
II is the same potential but with the running a s (r) given by Eq. ffTU]) . In this case, we 
have slightly changed the parameter M to preserve the value of V(oo,T c ) = 0.508 GeV 
that agrees with lattice estimate for the free quark-antiquark energy 0. Note that for the 
potential in Eq. ([9]) the Bargmann integral (j2J) is 

^(2-3—) (14) 



1 However, the difference in the parameter Mq causes the small difference of V^oo, T) for T > T c , see 
tables 03 El 
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The potentials III and IV are defined by the exact integral representation (jSJ) for V(r) 
and correspond to the constant and running a s , respectively. 

We display in Figs. [TJ |2]the binding energies of the IS J ftp and T mesons above the 
deconfinement temperature. The details of the calculation are presented in Tables [21 [3][E 
In these tables we present the constituent quark masses \iq for cc and 66, the binding 
energies E — V np (oo,T), the mean squared radii r = V< r 2 >, and the masses of the 
QQ mesons. We employ m c = 1.4 GeV and m;, = 4.8 GeV. As in the confinement region, 
the constituent masses /iq only slightly exceed the bare quark masses rriQ that reflect 
smallness of the kinetic energies of heavy quarks. We also mention that the account of 
the running a s (r) in the Coulomb potential produces a tiny effect as compared with the 
case of a constant a s = 0.35 both for the energies (compare lines I and II and III and IV 
in Tables El [3]) and 66 wave functions in Fig. [3j Both for the charmonium and for the 
bottomonium states the energy gap V(oo,T) — Eq(T) gets smaller and the mean square 
radius r gets larger as the temperature grows. We find no excited states IP and 2S 
states, although the unbound 2S 66 state appears to be very close to the threshold. 

At T = T c we obtain the weakly bound cc state. The potential II predicts small 
additional binding ~ 20 MeV as compared with the binding for potential I. At T = 1.3 T c 
the difference of V(oo, T) — E (T) calculated for different potentials comprises only a few 
MeV. The melting temperature for the case IV is ~ 1.3 T c . The charmonium masses lie 
in the interval 3.1-3.3 GeV, that agrees with the results of [T2]. Note that immediately 
above T c the mass of the cc state is about 0.2 GeV higher than that of J/ip. 

As expected, the T state is much more bound and remains intact up to the larger 
temperatures, T ~ 2.3 T c (all the details of calculation can be inferred from Table |3]). 
The masses of the L = bottomonium lie in the interval 9.7-9.8 GeV, about 0.2-0.3 
GeV higher than 9.460 GeV, the mass of T(IS') at T — 0. At T = T c the 66 separation 
tq is 0.25 fm, which is compatible with r$ = 0.28 fm at T = 0. At the melting point 
ro — > oo. Note that the IS bottomonium undergoes very little modification till T ~ 2T C . 
The results agree with those found previously for a constant a s = 0.35 |5J. The melting 
temperatures for the J/ip and T are shown in Table HI The results for the T are in 
agreement with the lattice study of Ref. [16j. Note that in our calculations we neglect 
the spin-spin force, therefore the J/ip and rj c mesons appear degenerate, as well as the T 
and rjb. This degeneracy is expected to be removed by a short range spin-spin interaction, 
whose effect, at T = 0, is often treated perturbatively assuming a contact interaction. 

4.2 Dissociation of J/ip and T in collision with gluons 

Heavy quark bound states are important probes of the dynamics in the QGP. Charmonium 
suppression has been observed at a variety of energies at SPS [T7] and RHIC [TS|. While 
the melting of bound states certainly reduces quarkonium production, the converse is 
not true: different, even competing effects make it difficult to interpret charmonium 
suppression patterns. It has been noted that such effects should be less significant for 
bottomonium [j~9] . 

2 In tables [5J [3] the results for the potential I that were previously reported in [S] are quoted for 
comparison. 
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Previous treatment of the dissociation of heavy quarkonium by the absorption of a 
gluon was carried out [20] using the operator product expansion in the the large N c limit 
and hydrogen states to evaluate the transition matrix elements. In heavy quarkonia of 
interest, the radial dependence of the quark-antiquark potential and the corresponding 
wave functions differs from the Coulomb potential. The calculation of the dissociation 
cross section for a color El transition can be well described by the potential model, follow- 
ing the results of Akhiezer and Pomeranchuk [21 J and Blatt and Weisskopf [22] obtained 
for the photo-disintegration of a deuteron. At low energies, the dominant dissociation 
cross section is the El color-electric dipole transition for which the quark final state will 
be the continuum (QQ) IP state. We calculate the cross section for the quarkonium 
dissociation after a gluon impact similarly the calculation of the deuteron disintegration 
via the photon absorption as has been done in Refs. 



An initial bound (QQ)\s state with a binding energy e(T) = V(oo, T) — E (T) relative 
to the threshold in the color-singlet potential □ absorbs a El gluon of energy u>, and is 
excited to the color-octet final state (QQ)ip with the energy 

— = u -e(T). (15) 



The corresponding cross section reads 

^QQ)» = ^^^T^I / "j/>(>')"i>'(>'. T)nlr I . (16) 



gQ K -f Kg 



where uis(r,T) is the wave function of the QQ bound state at the temperature T nor- 
malized as 

oo 

r \u ls (r,T)\ 2 dr = 1, (17) 



a 



gQ — a s/Q, Uip(r) is the free P-wave QQ wave function 

sin(Ar) 

Uip(r) = — cos(kr), 

kr 



k is the momentum of outgoing quarks in Eq. (115]) . and 

fc 2 = 2 yUQ e(T). (19) 

The results are shown in Figs HI [5] They are in a qualitative agreement with those of Ref. 
[23] where the QQ potential was identified with the color singlet heavy quark free energy 
above T c taken from quenched lattice QCD simulations. 



4.3 QQQ baryons at T > T c 

The three-quark potential is given by 

Vqqq — 



\Y1 V QQ^ T )i ( 2 °) 

i<j 



e(T) > below the melting point. 
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where | is the color factor. We solve the three-quark Schrodinger equation by the hyper- 
spherical harmonics method. Using the three-body Jacobi coordinates 




(ri - r 2 ) 



2 (r l + r 2 

f^Q I o ^3 



2 v . .. V 3 r v V 2 

the wave function i/i(p, A) in the hypercentral approximation is written as 

1 u{R,T) 



3 i?5/2 



where the hyperradius 



*(p,A,r) 



D 2 2 i \2 A'-Q / 2 2 2 \ 

# = p + A = — (r 12 + r 23 + r 31 J 



(21) 



(22) 



(23) 



is invariant under quark permutations, and 

p = Rsm9, \ = R cos 9, 0<6<n/2. (24) 

Averaging the three-quark potential fl20|) over the six-dimensional sphere, one obtains the 
one-dimensional Schrodinger equation for the reduced function u(R,T): 



d 2 u(R,T) 
d~R 2 



+ 2 



Eo-^- 2 -^(Vc(R,T) + V(R,T)) 



u(R,T) = 0, (25) 



where 



tt/2 

V C (R,T) = -\a s [ exp(-m d (T)R) 



R 



(26) 



V(R,T) = V(oo,T)- 



tj{T)B 



tt/2 



T 



V(oo,T) being given by Eq. ([7]), and 

M R sin 9 



R 



K 1 (R)R-—e- H (l+R))dn e , 

Mq 
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dn 6 = — sin 2 9 cos 2 9 d9. 



(27) 



The temperature-dependent mass of the colorless QQQ states is defined as 

3m 2 Q 3 
M QQQ = o + o^Q + Eo(m Q ,n Q ), 

where \iq are now defined from the extremum condition imposed on Mqqq in (129 

9M QQQ = Q 
9pq 

Note that the average interquark distances are 



(28) 



(29) 



(30) 



< rfi > 



< R 2 > 



(31) 
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see Eq. (I2"3"j) . The bound QQQ state exists if Eq^itiq, hq) < Vqqq(oo,T), where 

V QQQ (oo,T) = ?V(oo,T). (32) 

In our three-quark calculations we use the potential I. For this potential there is no bound 
^ccc states P. However, in all our calculations the Vt ccc was found to lie almost at threshold. 
The VL bbb survives up to T ~ 1.8 T c , see Tableland Fig. 00. 

5 Conclusions 

In conclusion, we have calculated binding energies and melting temperatures for the lowest 
eigenstates in the cc, bb, and bbb systems. The color electric forces due to the nonconfining 
correlator Df survive in the deconfined phase and they can support bound states at 
T > T c . For what concerns the charm states, we find that J/ip survive up to T ~ 1.3 T c , 
and there is no bound Q c state at T > T c . On the other hand, the bb and bbb states 
survives up to higher temperature, T ~ 2.6 T c and T ~ 1.8 T c , respectively. This suggests 
that the systems are strongly interacting above T c . 

This work was supported in part by RFBR Grant 09-02-00629. 
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Table 1: Parameters of the non-perturbative quark- ant iquark potentials described in the 
text. B and M are in units GeV 3 and GeV, respectively. The potentials I-IV correspond 
to the value of V(oo,T c ) = 0.508 GeV 





V{r) 


a l 


az s (r) 


B 


M 


I 


Eq. ® 


0.6 


0.35 


0.583 


0.9 


II 


Eq. © 


0.563 


Eq.dTO]) 


0.494 


0.813 


III 


Eq. OBJ 


0.6 


0.35 


0.583 


0.9 


IV 


Eq. (ED 


0.563 


Eq.flTDD 


0.494 


0.813 



Table 2: Details of the calculation of the cc states as a function of the temperature above 
the deconfmement region. V(oo,T), /i c , E — V np (oo), and M cZ are given in units GeV, r 
in units GeV -1 . m c = 1.4 GeV 



T/T c 


Potential 


V(oo,T) 




E - V np (oo) 


^0 


M c - C 


1 


I 


0.508 


1.451 


- 0.019 


7.53 


3.291 




II 


0.509 


1.469 


- 0.040 


6.07 


3.271 




III 


0.508 


1.454 


- 0.022 


7.24 


3.288 




IV 


0.509 


1.473 


- 0.048 


5.75 


3.264 


1.3 


I 


0.381 


1.419 


+ 0.006 


10.50 


3.186 




II 


0.372 


1.425 


+ 0.001 


9.75 


3.173 




III 


0.381 


1.424 


+ 0.001 


9.91 


3.183 




IV 


0.372 


1.434 


- 0.008 


8.72 


3.165 


1.6 


IV 


0.262 


1.416 


+ 0.007 


10.74 


3.069 
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Table 3: Details of the calculation of the bb states as a function of the temperature above 
the deconfmement region. The notations are the same as in Table El m& = 4.8 GeV 



T IT 
J- / J-c 


JrOlcIlCldl 


V [OO, 1 ) 




TP T/ f/-v~A 


TO 


M - 
M bb 


1 


T 

I 


0.508 


4.984 


- 0.300 


1.27 


9.815 




TT 
11 


U.OUo 


4.90O 


- U.olo 


1 QO 

l.oz 


9. CJo 




TTT 
111 


0.508 


4.985 


- 0.308 


1.26 


9.807 




IV 


0.509 


4.953 


- 0.326 


1.31 


9.786 


1.3 


I 


0.381 


4.950 


-0.183 


1.55 


9.802 




TT 
11 


n Q7o 
U.o ( 2 


A no K 


- U.lo ( 


l.Dl 


y. loo 




TTT 
111 


0.381 


A n r o 

4.953 


n onn 

- 0.200 


1.51 


9.785 




IV 


0.372 


4.928 


- 0.211 


1.55 


9.764 


1.6 


T 

I 


0.275 


4.915 


-0.095 


2.06 


9.783 




TT 
11 


U.ZOZ 


4.o9o 


- U.U9o 


Z. 14 


9. M 2 




TTT 
111 


n o*r r 

0.275 


A H0 1 

4.921 


A 1 1 n 

- 0.119 


1 n 1 

1.91 


9.759 




IV 


0.262 


4.902 


A "1 O/' 

- 0.126 


"i n /i 

1.94 


9.739 


2.0 


I 


0.162 


4.863 


-0.021 


4.25 


9.742 




TT 
11 


U.140 


4.001 


- U.U1 ( 


4.DD 


n von 

y. r /y 




TTT 
111 


n 1 

0.162 


4.878 


n n /i ^ 

- 0.046 


o no 

3.03 


n T1 T 

9.717 




1 V 


U.140 


4. ODD 


n n/i s 
- U.U4o 


O.UO 


y.oyo 


2.2 


T 
1 


n i i r 

0.115 


/I o o o 

4.832 


n nno 

- 0.003 


7.62 


n "r i o 

9.712 




TT 
11 


U.uy ( 


A fiOQ 
4.0ZO 


n nm 
- U.UUl 


o.Ol 


y.oyo 




TTT 
111 


0.115 
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4.847 


-0.023 


4.34 
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I 


0.093 


4.818 


+ 0.001 


9.52 


9.694 




II 


0.075 


4.813 


+0.002 


10.17 


9.677 


2.4 


III 


0.073 


4.831 


-0.007 


6.80 


9.665 




IV 


0.054 


4.827 


-0.007 


6.91 


9.648 


2.6 


III 


0.034 


4.812 


+0.001 
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IV 


0.016 


4.809 


+0.001 


10.18 
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Table 4: Dissociation temperatures (in units of T c ) for cc and bb states 
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II 


III 


IV 




1.24 


1.29 


1.29 


1.46 


T 


2.27 


2.29 


2.57 


2.57 
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Table 5: Details of the calculation of bbb baryon above the deconfinement region for 

the potential I of Table [TJ Vqqq(oo,T) is defined by Eq. (132)) . \/~B? is related to the 
interquark distances by Eq. (13T]) . Dimensions are the same as in Table [2j m;, = 4.8 GeV 



T 

T r 


Vqqq(oo,T) 




Eq- V(oo,T) 


\[~R 2 


M bbb 


1 


0.763 


4.943 


- 0.279 


3.73 


14.890 


1.3 


0.571 


4.906 


-0.135 


4.80 


14.840 


1.6 


0.414 


4.865 


-0.037 


7.17 


14.777 


1.8 


0.324 


4.835 


+ 0.001 


10.00 


14.724 
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Figure 1: V(oo,T) (thick curve) and E Q (T) (thin curve) are plotted against T/T c for the 
J/ip (the potential IV). V(oo,T) and E (T) are given in units of GeV. 
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;ure 2: The same as in Fig. [I] for the ground T state 
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Figure 3: The bb wave function u(r,T) normalized to unity (r in units of GeV -1 ). Solid 
and dashed curves show the u(r,T c ) for the potential I and II, respectively, dot-dashed 
and dot curves show the u(r, 2T C ) for the potentials I and II. 
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Figure 4: Cross sections a for J/ip dissociation in the potential model vs gluon energy ou. 
The curves are labelled by the ratios T/T c . a and u are given in units GeV~ 2 and GeV, 
respectively. 
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Figure 5: Cross sections for bb dissociation in the potential model. The notations are the 
same as in Fig. H] 
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Figure 6: V(oo,T) (thick curve) and E (T) (thin curve) in units of GeV for the Vt hbb (the 
potential I) 
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